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Abstract 

We report on our results about hadronic tt production at NLO QCD including t, t 
spin effects, especially on tt spin correlations. 

1. Introduction 

Top quarks, once they are produced in sufficiently large numbers, are a sensitive 
probe of the fundamental interactions at high energies. On the theoretical side this 
requires precise predictions, especially within the Standard Model (SM). As far 
as tt production at the Tevatron and LHC is concerned, spin-averaged (differen- 
tial) cross section have been computed at next-to-leading order (NLO) QCD [1,2] 
including resummations [3,4]. Observables involving the spin of the top quark 
can also be calculated perturbatively, especially within QCD. It is expected that 
such quantities will play an important role in exploring the interactions that are 
involved in top quark production and decay. Within the SM, the QCD-induced 
correlations between t and t spins are large and can be studied at both the Tevatron 
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and the LHC, for instance by means of double differential angular distributions of 
tt decay products. Results for these distributions at NLO QCD [5-8] are reviewed 
below. 

2. Theoretical Framework 

We consider the following processes at hadron colliders 

r /+/'- + x 
pp/pp^tt+x^ ! + x (i) 

{ hit + x 

where I (/') = e,/J,x, and j t (jj) denotes the jet originating from non-leptonic t (F) 
decay. At NLO QCD the following parton reactions contribute to the above pro- 
cesses: 

sg,qq bb+Af, 
gg,qq bb+Af+g, 

g + q (q)JUbb + 4f + q(q), (2) 

where f = q 1 £,V£. The calculation of these cross sections at NLO QCD simplifies 
in the leading pole approximation (LPA), which is justified because T t /m t < 1%. 
Within the LPA, the radiative corrections can be classified into factorizable and 
non-factorizable contributions. At NLO these non-factorizable corrections do not 
contribute to the double differential distributions [8], which we will discuss be- 
low. Therefore we do not consider them here. Considering only factorizable 
corrections in the on-shell approximation the squared matrix element \M | 2 of the 
respective parton reaction is of the form 

\M | 2 oc Tr [p#p] = pa'atf«a',pp' P(3p'- (3) 

Here R denotes the density matrix that describes the production of on-shell tt 
pairs in a specific spin configuration, and p, p are the density matrices describing 
the decay of polarized t and t quarks, respectively, into specific final states. The 
subscripts in © denote the t, t spin indices. The spin-averaged production density 
matrices yield the NLO cross sections for tt being produced by qq, gg, gq, and gq 
fusion [1,2]. 
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To obtain a full NLO QCD analysis of <UJ), we must consider also the NLO QCD 
corrections to the matrix element of the main SM decay modes of the (anti)top 
quark in a given spin state, i.e. the semileptonic modes t — > b£ + \£, b£ + Vgg (£ = 
e,fi,%), and the non-leptonic decays t — > bqq' ', bqq'g where qq' = ud,cs for the 
dominant channels. For the computation of the double angular distributions ©, 

the matrix elements of the 2-particle inclusive parton reactions i — > a + b+X are 
required. Here a,b denote a lepton or a jet. In the LPA this involves the 1 -particle 
inclusive t decay density matrix 2p^ a = r^^ll + KaT ■ qi) a ' a , where qi is the 
direction of flight in the t rest frame and is the partial width of the respective 
decay channel. An analogous formula holds for t decay. The factor K a is the t 
spin analysing power of particle/jet a. Its value is crucial for the experimental 
determination of top spin effects, in particular of it spin correlations. For the 
standard V — A charged-current interactions these coefficients are known to order 
a s for semileptonic [9] and non-leptonic [10] modes. The charged lepton is a 
perfect analyser of the top quark spin, which is due to the fact that K| = 1 — 
0.0150$. In the case of hadronic top quark decays, the spin analysing power of 
jets can be defined, for example 

Kh = -0.408 x (1 -0.340o 5 ) = -0.393, (4) 
Kj = +0.510 x (1 -0.654O,) = +0.474. (5) 

Here is the analysing power of the b jet and K 7 refers to the least energetic 
non-^-quark jet defined by the Durham algorithm. Obviously the spin analysing 
power is decreased if one uses the hadronic final states to analyse the spins of t 
and/or t. However, this is (over)compensated by the gain in statistics and by the 
efficiency with which the t (J) rest frames can be reconstructed. 
With the above building blocks, we can discuss the following double angular dis- 
tributions 1 

1 do 1 / \ 

—a ^a 5" = 7 1-Ccos0icos0 2 , (6) 

Gacos0iacos02 4 \ / 

where C is a measure of the it correlations; 0i (62) is the angle between the direc- 
tion of flight of particle/jet a\ [ai) in the t (F) rest frame with respect to reference 
directions a (b), which will be specified below. For the factorizable corrections 
the exact formula C = K ai K fl2 D holds [7]. Here D is the it double spin asymmetry 

_ iV(TT) +iV(|l) -A^(TI) -iV(lT) 

jv(TT)+jv(u)+jv(ti)+;vatr 



'QCD-generated absorptive parts in the parton scattering amplitudes induce a small t and t 
polarization, which to order oc^ is normal to the qq,gg — > it scattering planes [12, 13]. 
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where iV(TT) denotes the number of ti pairs with t (t) spin parallel to the reference 
axis a (b), etc. Thus a and b can be identified with the quantization axes of the 
t and t spins, respectively, and D directly reflects the strength of the correlation 
between the t and t spins for the chosen axes. 

For ti production at the Tevatron it is well known that the so-called off-diagonal 
basis [11], which is defined by the requirement that 6(tl) = 6(IT) = f° r me 
process qq — > ti at tree level, yields a large coefficient D. It has been shown in [7] 
that the beam basis, where a and b are identified with the hadronic beam axis, is 
practically as good as the off-diagonal basis. A further possibility is the helicity 
basis, which is a good choice for the LHC. 

3. Predictions for the Tevatron and the LHC 

We now discuss the spin correlation coefficients C of the distributions ©. It 
should be noted that beyond LO QCD, it is important to construct infrared and 
collinear safe observables at parton level. In the case at hand it boils down to the 
question of the frame in which the reference directions a and b are to be defined. 
It has been shown that, apart from the t and i rest frames, the ti zero momentum 
frame (ZMF) is the appropriate frame for defining collinear safe spin-momentum 
observables. The off-diagonal, beam, and helicity bases are defined in the ti ZMF. 
Details can be found in Ref. [8]. 

In Table[l]we list our predictions for C in © at the Tevatron and LHC. The results 
are obtained using the CTEQ6L (LO) and CTEQ6.1M (NLO) parton distribution 
functions (PDF) [14]. Numbers are given for the dilepton (L-L), lepton+jet (L-J) 
and all-hadronic (J-J) decay modes of the ti pair, in the latter two cases the least 
energetic non-Z?-quark jet (defined by the Durham cluster algorithm) was used as 
spin analyser. One notices that for the Tevatron the spin correlations are largest in 
the beam and off-diagonal bases, and the QCD corrections reduce the LO results 
for the coefficients C by about 10% to 30%. For the LHC the QCD corrections are 
small (< 10%). These results are obtained with ^ = /jr = /up = m t = 175 GeV. At 
the Tevatron, a variation of the scale /j between m t /2 and 2m t changes the results 
at jU = m t by ~ ± (5-10)%, while at the LHC the change of Ch e i is less than a per 
cent. 

In Table 12 we compare the NLO results for the spin correlation coefficients eval- 
uated for the CTEQ6.1M, MRST2003 [15] and GRV [16] PDFs. It is easy to see 
that the results with the recent CTEQ6.1M and MRST2003 PDF agree at the per 
cent level (this is not the case for previous versions of the CTEQ and MRST PDF), 
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Table 1 : LO andNLO results for the spin correlation coefficients C of the distribu- 
tions $@j)for the Tevatron at y/s = 1 .96 TeV and for LHC at y/s =14 TeV. The PDF 
CTEQ6L (LO) and CTEQ6.1M (NLO) were used, and /j F =/j R = m t = \15GeV. 







L-L 


L-J 


J-J 


Tevatron 


Chel 


LO 
NLO 


-0.471 - 
-0.352 - 


-0.240 
-0.168 


-0.123 
-0.080 


Cbeam 


LO 
NLO 


0.928 
0.777 


0.474 
0.370 


0.242 
0.176 


C ff 


LO 
NLO 


0.937 
0.782 


0.478 
0.372 


0.244 
0.177 


LHC 


Chel 


LO 
NLO 


0.319 
0.326 


0.163 
0.158 


0.083 
0.076 



while the GRV98 PDF gives significantly different results at the Tevatron. This 
shows that the spin correlations are very sensitive to the relative quark and gluon 
contents of the proton [7]. Future measurements of © may offer the possibility 
to further constrain the quark and gluon contents of the proton. 
Before closing this section, we summarize how an experimental measurement of 
the distributions © that matches our predictions should proceed: 1) Reconstruct 
the top and antitop 4-momenta in the laboratory frame (= cm. frame of the 
colliding hadrons). 2) Perform a rotation-free boost from the laboratory frame 
to the tt ZMF Compute a and b in that frame. 3) Perform rotation-free boosts 
from the tt ZMF to the t and t quark rest frames. Compute the direction qi (%) 
of the t (F) decay product a\ (a.2) in the t it) rest frame. Finally, compute cosGi = 
a ■ qi, cos 82 = b • q2. Here one should notice that in this prescription the t and t 
rest frames are obtained by first boosting into the tt ZMF. If this step is left out, 
and the t and t rest frames are constructed by directly boosting from the lab frame, 
a Wigner rotation has to be taken into account. 
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Table 2: Spin correlation coefficients at NLO for different PDF s for the Tevatron 
( upper part) and the LHC (lower part) for dilepton final states. 



Tevatron 




CTEQ6.1M MRST2003 


GRV98 


Chel 


-0.352 -0.352 


-0.313 


Cbeam 


0.777 0.777 


0.732 


C ff 


0.782 0.782 


0.736 


LHC 


Chel 


0.326 0.327 


0.339 



4. Conclusion 

We have computed at NLO QCD the tt spin correlations in hadronic top pro- 
duction, which are large effects within the SM. Our present results are obtained 
without imposing kinematic cuts. Such cuts will in general distort the distribu- 
tions, i.e. C will in general depend on the angles 0i and 02- One strategy is to 
correct for these distortions by Monte Carlo methods before extracting the spin 
correlation coefficient and comparing it with theoretical predictions. A future aim 
is to directly include the cuts in an NLO event generator to be constructed with 
our NLO results for all relevant 2^6 and 2^7 processes. 
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